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We study the diffusion of a Brownian particle quadratically coupled to a thermally fluctuating field. In the 
weak coupling limit, a path-integral formulation allows to compute the effective diffusion coefficient in the 
cases of an active particle, that tends to suppress the field fluctuations, and of a passive particle, that only 
undergoes the field fluctuations. We show that the behavior is similar to what was previously found for a lin- 
ear coupling: an active particle is always slowed down, whereas a passive particle is slowed down in a slow 
field and accelerated in a fast field. Numerical simulations show a good agreement with the analytical cal- 
culations. The examples of a membrane protein coupled to the curvature or composition of the membrane 
are discussed, with focus on the room for anomalous diffusion. 



I. INTRODUCTION 

Diffusion of an object interacting with its fluctuating 
environment has recently received a lot of attention: ex- 
perimental studies have investigated the cases of colloidal 
beads diffusing along lipid bilayer tubes or through an actin 
network |1|, insulin granules diffusing in /3-cells |2|, or di- 
electric colloids subject to random optical forces generated 
by multiply scattered light [3] . A variety of behaviors have 
been observed, most of the time including anomalous dif- 
fusion, either in the mean squared displacement or in the 
probability distribution function. These observations call 
for a general theoretical framework able to describe diffu- 
sion in a complex environment. The considered systems 
can be cast into two classes: if the object affects its environ- 
ment, as in the first and second examples, it is called active , 
if it does not, as in the third example, it is called passive . 

Among all the investigated systems, membrane proteins 
have concentrated much effort. One of the first theoreti- 
cal studies is due to Saffman and Delbriick |4|, who have 
computed the hydrodynamic drag felt by a protein moving 
at constant velocity in a membrane; using the Einstein re- 
lation (5], it allows to determine the diffusion coefficient. 
However, this calculation only gives a weak logarithmic de- 
pendence of the diffusion coefficient on the protein size that 
was contradicted later by accurate experiments (6| . Explain- 
ing the observed diffusion coefficient of membrane proteins 
has thus remained a major theoretical challenge. The nu- 
merous analytical and numerical investigations of the effec- 
tive diffusion coefficient aim at taking into account two new 
effects: the geometrical penalty of free diffusion on a ruffled 
surface 1 7 -12 1 and the interaction of the protein with a local 
parameter, most of the time the membrane curvature [TlfTOl- 
[T2l but also its height (T3|. These works have shown that 
the protein diffusion coefficient is reduced, whereas it is in- 
creased if the protein action on the membrane is neglected 
im [121 . This result is derived when the membrane equili- 
brates much faster than the protein moves; this is called the 
adiabatic limit and it is the biologically relevant one 1 12 1 . 

In previous studies fTiHlTl , we did a move towards a more 
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FIG. 1: (Color online) Particle trajectory in a fluctuating field; it is 
quadratically coupled to the field and sees the effective potential 
|(/.(3/,f)2. 



general model, that aims at describing not only membrane 
proteins but also every kind of objects moving in a fluctuat- 
ing environment, such as the ones presented sooner; a very 
generic picture is given in FIG[T] In this model, a particle 
moving in a space of arbitrary dimension is linearly cou- 
pled to a scalar Gaussian field. We first computed the drag 
force felt by the object when it is pulled at constant velocity 
[14l[15l. Naively, the Einstein relation would allow to de- 
duce the diffusion coefficient but this is not the case since 
the drag is computed at constant velocity whereas the Ein- 
stein relation needs the drag at constant force . To clarify this 
issue, we computed the diffusion coefficient in the adiabatic 
limit, where the field is much faster than the particle, and 
showed that it can indeed be inferred from the drag coeffi- 
cient computed at constant velocity [16, 17J. In this limit, 
we showed that an active particle is always slowed down by 
its coupling to the field, whereas a passive particle is always 
accelerated: this corresponds to what was observed for the 
motion of membrane proteins |TT] [12] . However, in a more 
general theory one has to consider the non-adiabatic limit, 
where we showed that the previous statement breaks down: 
in a slow field a passive particle is slowed down too [17] . 

This first generalization needs to be further extended in 
order to be able to model more couplings between the ob- 
ject and the environment. For instance, when a protein is 
coupled to membrane curvature, there is a superposition of 
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two effects: the protein may impose a spontaneous curva- 
ture to the membrane, and, if it is stiffer than the membrane, 
it can also reduce the membrane fluctuations. Whereas the 
first effect is well described by a linear interaction between 
the protein and the membrane, the second one needs a 
quadratic interaction. What was proved for the linear in- 
teraction needs to be investigated for the quadratic inter- 
action: this article fills this gap. The drag force at constant 
velocity has already been computed for the quadratic inter- 
action [18] . 

In this article, we introduce the general model in section 
[njand explain how it describes a protein coupled to mem- 
brane curvature or composition. Analytical calculations are 



performed in section III Firstly, we derive the exact path- 



integral representation of the particle motion. Secondly, we 
show that in the weak coupling limit it reduces to an effec- 
tive dynamics for the particle, containing a memory term. 
Thirdly, we address the adiabatic limit, derive a Markovian 
dynamics and compute the effective diffusion coefficient; it 
appears that it can be deduced from the drag computed at 
constant velocity. Fourthly, we compute the effective diffu- 
sion coefficient for a weak coupling but outside the adia- 
batic limit and show that an active particle is always slowed 
down whereas a passive particle is slowed down in a slow 
field and accelerated in a fast field. These results for a 
quadratic coupling, albeit very similar to those for a linear 
coupling, are new; moreover, to get to the effective diffusion 
equation, the quadratic coupling requires new techniques 
that are detailed here. Analytical results are compared to 
numerical simulations on a simple unidimensional system 
The possibility for anomalous diffusion is dis- 



m section 



IV 



cussed in sectionM 



II. MODEL 



A. General model 



Our general model is relatively simple but contains sev- 
eral adjustable parameters that allow to use it to model 
many different systems. Some of these ones are given as ex- 
amples at the end of this section. 

This model is close to the one presented in (17), with the 
exception that the coupling with the field is here quadratic 
in the field, making the interaction "fluctuations induced". 
This quadratic coupling has already been introduced in (18l , 
albeit a in less general formulation, where the particle was 
pulled at constant velocity; here the particle is submitted to 
a thermal noise. 

We consider a particle located at x in a d-dimensional 
space, with a quadratic coupling to a Gaussian field </)(y). 
The energy of the particle-field system is given by 

H[(p, x]^^jcl)[y] lA(p] (y)dy + ^ [Kcp] [xf. (1) 



field, and the operator A gives the shape of the field we con- 
sider. The second term is the quadratic coupling between 
the field and the particle, the particle shape is defined by 
the operator K, and h is the coupling constant. The pure 
Casimir interaction where the particle suppress the field 
fluctuations is obtained in the limit h ^ oo. The opera- 
tors A and K are translation invariant and isotropic (cf. ap- 
pendix|A}. 

We start with an overdamped Langevin equation for the 
whole system, the equilibrium of which is thus given by the 
Gibbs-Boltzmann statistics. The particle position evolves as 
follows, 

SH 



Xit) = -Kjc—lcpiy, t),X[t)] + sfK^T]{t], (2) 

ox 



hXy 



V[(Kcl)f]{x{t),t] + ^r](t]. 



(3) 



Kx is the mobility of the particle and i/ ( t) is a Gaussian white 
noise with correlator 



{v{t)Tj[tr)^2TSit-t'n. 



(4) 



The superscript "T" denotes the transposition and 1 is the 
identity matrix. The field evolution is 



+ h{K(p]{x{t],t] {RK5] {y-x{t]]] 
+ ^[/R£)^ {y,t). 



(6) 



is the "mobility" (or evolution rate) of the field, i? is a 
translation-invariant isotropic operator defining the field 
dynamics, and ^(y, t) is a functional Gaussian white noise 
with correlator 



(ay, my', t')) ^2T5[y-y')S{t- t']. 



(7) 



To take into account the case of a passive particle, i.e. a par- 
ticle that does not affect the field, we introduce a feedback 
parameter C in j6| that becomes 

(/)(y, f) = -K0[(i?A0)(y, t) 

+ (h{K(l)){x[t), f) [RK6) [y-xm] 

+ ^[VRi]{y,t). (8) 

The previous equation is recovered with ( - I whereas a 
completely passive particle corresponds to ^ = 0. Except for 
( - I, the dynamics does not satisfy detailed balance, and 
the system is out of equilibrium. A trajectory of the particle- 
field system is given as an example in FIG.[T] 

In this paper, we will focus on the effective or late -time 
diffusion coefficient, defined by 



The notations used for functional operators are defined in 
appendix|A| The first term is the quadratic energy of a free 



(9) 
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Of course, this number is only defined if tfie diffusion is nor- 
mal, and we will restrict ourselves to this case here. How- 
ever, it can give a hint on the kind of diffusion: Des - 
indicates subdiffusion and Deft - oo signals superdiffusion. 
Without coupling, the bare diffusion coefficient is 

(10) 



B. Examples 

Our model applies to many systems where the presence 
of an inclusion penalizes the fluctuations of an order pa- 
rameter in its environment. In our examples the inclusion is 
a membrane protein that can be coupled to the membrane 
curvature or composition. The dimension of the space ac- 
cessible to the protein is d-2. 

The case of a protein affecting the curvature has been ex- 
tensively discussed |7-12 19 20 1. First, we define the free 
field quantities: the membrane height corresponds to our 
field <p{y, t) and its free field energy is given by 

Hm^\^ \Km{^^ip{y)f + a{V(p{y)f\dy, (11) 

where is the membrane bending modulus and a its sur- 
face tension. In our formalism, the corresponding A opera- 
tor is given in Fourier space by 



(12) 



where m - y/alK^a is the inverse of the correlation length 
/. The protein tries to impose its spontaneous curvature, 
giving rise to the coupling 



H;r,t{<P,X] = ^ [V2(/)(x)-2Cp]'- ^ [V^<p{x)f, 



(13) 



where Cp is the spontaneous protein curvature and TCp its 
bending rigidity. Although in most of the studies the inter- 
action has a non-trivial shape, we need to consider it point- 
like here. The size of the protein can be introduced as an 
effective cut-off at the end of the computation. This inter- 
action term can be decomposed into a linear contribution, 
that tends to deform the membrane, plus a quadratic con- 
tribution, that tends to suppress its fluctuations. The case of 
a linear interaction has been treated in \ Y7\, and here we are 
only interested in the quadratic term. Physically, this corre- 



sponds to a tough, but flat protein: Atc = 
and the interaction is defined by 

h = Atc, k[k) = 



fp-'fm>0, Cp = 0, 



(14) 



Finally, the membrane dynamics is given by the Oseen hy- 
drodynamic tensor |21 1, 



R{k)^ 



47?|7c| 



(15) 



where 77 is the viscosity of the surrounding fluid. The effect 
of projection due to the diffusion on a ruffled surface I7HT21 
is neglected here. 



Lipid membranes are made of several kinds of lipids, and 
proteins can also be coupled to the composition field (p\22\, 
that fluctuates on growing length scales as the system is 
tuned close to the miscibility transition |23|. The compo- 
sition field is ruled by the follov«ng operators 122] 



A(fc) = Jt^-hm^, 
K(k) = 1, 
R(k) = k^. 



(16) 
(17) 
(18) 



The operator R ensures that the average composition is con- 
served. Again, m is the inverse correlation length, and goes 
to zero close to the miscibility transition as 



m ■■ 



~ [l-TITcf, 



(19) 



where Tc is the miscibility temperature and v = 1 is the ex- 
ponent observed experimentally on Giant Unilamellar Vesi- 
cles |23|. The interaction between two proteins has been 
calculated exactly at the critical point in 1 24| . 

Proteins can be coupled to other membrane fields such as 
the local membrane thickness 1 13|. These parameters can 
be advected, or not, by the flow created by the protein mo- 
tion. This effect, discussed in (25l|26j, is neglected here. 



III. ANALYTICAL CALCULATIONS 
A. Path-integral formalism and effective action 

The evolution of the particle-field system, that is com- 
pletely defined by j3]|4][7][8}, can be cast in a path-integral 
formulation. Following the procedure used in | IT, 27| with 
the Ito convention, the partition function for the whole sys- 
tem is of the form 

Z = J exp So Ix, p] - [(/), y/] - Sint [x, p,(p,f]^ 

X [dx]ldp][d(l)]ldy/]. (20) 

pit] and y/(.y, t) are the response fields |27| associated re- 
spectively with xit) and (p{y, t), and the action contains 
three terms. So [x, p] is the action of the pure Brovmian mo- 
tion (i.e. without interaction with the field) and reads 1 17 

So[x,p] = -iJ p{t)-xit)dt + Dj,j pitfdt. (21) 

The free field action is of the same form: 

S^[(/),V] = -i j V[y,t)[<p{y, t) + K^{RAcp){y, t)]dydt 

+ Tk^ jy/{y, t){Ry/){y, t)dydt. 



(22) 



Finally, the interaction action couples the particle and the 
field: 

Sintlx,p,(l),yf]^-i^ j p{t]-V{[K(l}f){x{t],t)dt 

-i(hK(i, j {K(l)){x{t),t)[KRy/){x{t),t)dt. 

(23) 
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To get an effective action S?^ for the particle, we integrate 
over ((/), y/) : 

exp[-S^lx,p]]^ 

j exp(-S^[(/),i/A]-Sint[A:,p,0,i//])[d(^][di//]. (24) 

This integral is quadratic and can be computed. We need 
the following formula: for a random Gaussian vector X and 
a matrix A, 



(exp(x'^^X))^ = det(l-2^(XX'^)) ^'^ 



- exp 



;tr(log[l-2^(XXT)]) 



(25) 



(26) 



We use it considering [cp, if/) as a Gaussian random variable 
weighted by Sg , the interaction action Sint [x, p, (p, if] playing 
the role of X^AX. Fourier transforming the fields in space 
allows us to write the interaction action 



Sint[x,p,if>,y/] = j [Ai>^{-k,-k', t, t')4>{k, t)4>{k', t' 



+ Asw{-k,-k', t, t')(f>{k, t)Tif/{k', t')] ^^^f drdf', (27) 



where 



(28) 



A^^(k, k', t, r') = — ^Pf^) " (*^+ k')K{k)Kik') 

xe-'*+'^''-^t«5(f-f'), 

A^^{k,k',t, t')^-\i:hK^K{k)K[k')R[k')e-'^^^^'^-''^'^8{t- t'). 

(29) 

We then need the free field correlation functions that are 
computed in appendbc|B] 

C^0(fc, k', t, t') = (0(fc, t)4){k', t')) (30) 

= -J—e-''^^'-''^^*^'^*\2nf8{k+k'), (31) 
A(fc) 

C^y,{k,k', t, t') = (0(fc, t)if/{k', t')) (32) 
= ie-H\t-t'\m)K(k)g^^_ t')[2nf8{k+k'). 

(33) 

0(f) is the Heavyside function, with 0(0) - since we use the 
Ito convention. We introduce the operator 



— AipipC(p(p + A^^Cxjjip, 
whose Fourier transform reads 



(34) 



M(fc, k', t, t') = he-'^''^'' ■)-^(t)-K^\t-t'\R(k )Mk 
Tk,, p{t)-(k+k') 



A(fc') 



+ (K^R{k']e{t'-t) 



(35) 



The formula 26 gives the effective interaction action as 



1 



S!^lx.p]^-ti[\og{l + 2M[x,p])]; 



(36) 



it allows us to write the partition function 

Z = j exp(-So[x,p]-Sfft[x,p])[dx][dp]. (37) 
This expression is exact and constitutes our first new result. 



The logarithm appearing in the interaction action 36 



makes it very difficult to handle; its signification is discussed 
later. To go further, we have to expand the logarithm; this 
can be done for a small coupling h«l. 



B. Effective particle dynamics 

In the weak interaction limit (/i « 1) a simpler expression 
can be obtained by expanding the effective interaction ac- 
tion pe) in powers of the coupling constant h. We start by 
developing the logarithm it contains in powers of M, that is 
proportional to h: 

log(l + 2M)^2M-2M^ + &{h^). 



(38) 



Now, the trace can be taken explicitly: 



tr(M) = / M(fc,-7c,f,f) — ^df = 0, (39) 
J {2n)<^ 

since 0(0) = with our convention. The fact that the first 
term, proportional to h, does not contribute is not surpris- 
ing: to feel the field, the particle must interact with it at least 
twice. For the second term, we need 

trfM^]^ f M[k, k', t, t')M{-k',-k, t' , t) ^^^f, drdf' 
^ ' J {2n)^<^ 

{k+k')-\x(t)-x{t')]-K^\R(k)Klk)+R{k')Mk')]\t-t'\ 



dk 



h^TK^j e'' 



(K^p{t)-[k+k')^^6{t-t') 

_ TK^p[t)-[k+k')p{t')-{k+k') 
4 A{k)A[k') 

X K(kfK(k'f ^dtdt'. 

'■2d 



(2n) 



(40) 



This expression has been obtained after several variable 
changes on the different terms. The action at the order 
is thus of the form 

= i/ p{t)-F{x{t)-x{t'),t- t']dtdt' 
Jt>t' 

+ ^ j pitfG{x{t)-x{t'),t- t')p{t')dtdt', 



(41) 



where 



F{x,t)^ii;h^TK^K^j [k + k') 



K{kfK[k'fR{k) 



A(fc') 



X e 



i{k+k')-lxit)-xit')]-K,p[R[k)A{k)+R{k')A{k')]\t-t'\ ^fcdfc' 



{271)^'' ' 

(42) 
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and 



G(x, t) = — ^ / {k+k'){k+k 



,^j K[kfK[k'f 



Jl.k+k')-lxlt)-x{t')]-K,p [R{k)Mk)+R{k')Mk')]\t- 1'\ dfcdfc' 



(43) 



In the truncated interaction action 41 , we recognize the 



action associated to the effective evolution equation (TT) 



+ f F(x(t)-x(t'],t-t']dt' + E.(x(t),t], (44) 

J—oo 

where the noise S(x, f) has a correlator 

(S(x(0, f)2(x(f'), f')^) = rG(jc(r) - jc(r'), r- f'). (45) 

This dynamics is of the same form as the one obtained 
in [171 for a linear coupling between the field and the par- 
ticle. The case of a linear coupling being simpler, the ef- 
fective dynamics is exact, whereas it is a weak coupling ex- 
pansion here. This dynamics reveals a two-times interac- 
tion with the field, at times t' and t. This is also the sense 



of the /i-expansion performed in the effective action 36 to 



get the effective dynamics: to restrict ourselves to two-times 
interactions. Expanding further the logarithm in this ac- 
tion would add more terms to the effective dynamics, corre- 
sponding to four-times interactions, six-times interactions, 
etc. 



C. Adiabatic limit 

In the adiabatic limit, where the field is much faster than 
the particle « k^), the field memory time reduces to 
zero and the effective dynamics becomes Markovian. We 
follow the procedure used in 1 16 1 and compute F and G to 
the first order in KxIk^. For F, this amounts to write 

F[x{t)-x{t'),t- t') = 

(f- f')jfc(r)-VF(0,r- t']+G[[t- t'f). (46) 

Integrating over t' , we get to the first order in Kj^Ik^ 



/ 

J-C 



F[x{t)-x{t'),t- f')df' = 

/ 



Kip 2d 

{k^ + k'^)K[k)^K{k'f dkdk' 
A(fc) A(fc') [^(fc) A(fc) + R{k')Mk')] (2;r)2rf 



x{t). 



(47) 



For G, we write 

G[xit)-x{t'),t- t') = G(0,r- t') + 0[t- t'), (48) 



and we integrate over f' to get, again to the first order in 

K^lKip, 



K Tk 
G[x{t)-x{t'),t-t')^8{t-t')- 



I 



K(p d 
{k^ + k'^)K[kfKik'f 



dkdk' 



A(fc) A(fc') [^(fc)A(fc) + Rik')K{k')] (.271)^'^ 



(49) 



We recognize the same integral in 47 and ]49) , so we intro- 
duce 



{k^ + k'^)KikfKik') 



1^2 



dkdk' 



2d J A{k)A{k')[R{k)A{k) + R{k')Aik')] (.271)2'^ 



, (50) 



that allows us to write the effective evolution equation for 
the particle, to the order in the coupling and KjcIk^ in 
the evolution rates. 



x{t) = -C— Ax(f) + Vl^rj[t] + J^Ar]'{t], 

Kif, V Kip 



(51) 



where ij' ( f ) is a Gaussian white noise with the same correla- 
tor as r]{t) j4}. 

We recognize the drag coefficient found in 1 18|. This re- 
sult is very similar to the case of a linear coupling 1 16|: in 
the adiabatic limit, the drag coefficient allows to derive a 
Markovian effective dynamics; we have shovm here that this 
property also holds for a quadratic interaction, and may 
thus be universal. 

With this dynamics, the effective diffusion coefficient is 
easy to compute: 



adia 
eff 



= D;,|l-(2(-l)-^A|. 



(52) 



This expression indicates that the diffusion coefficient is de- 
creased for an active particle and increased for a passive 
one, and corresponds to what was found for proteins cou- 
pled to membrane curvature |11 12 1. We show in the fol- 
lowing that this is only true in the adiabatic limit. 



D. Diffusion coefiicient outside the adiabatic limit 

The effective diffusion coefficient can also be computed 
outside the adiabatic limit; this was done for a linear cou- 
pling in 1 16 1 for the active case and in |17| for the general 
case. The path-integral method used in 1 17 1 can be general- 



ized and applied to the effective diffusion equation 44 . 
First, we put the functions F and G in the general form 



Fix, r) = i J e^^'®-^-"'®l^l/(Q)dQ, 
G{x, t)^j e'^'®-^-'^'®l*lg(Q)dQ, 



(53) 
(54) 
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where in our case 



Q= (fc,fc'), 



fiQ) 
giQ) 



■■ [R{1c)A{1c) + R{1c')A{k')] , 
h^DjcK^ K{kfK{k'fRik) 

■ — ^r(Q) 



{2n) 
2{2n) 



2d 



A(fc') 
jK[kfK[k'f 
A(fc)A(fc') 



(55) 
(56) 
(57) 

(58) 
(59) 



In this form, the functions F and G are very similar to those 
found in fTTl for a linear coupling to the field. The treat- 
ment that leads to the effective diffusion coefficient is thus 
the same; we briefly recall it but refer the reader to 1 17| for 
more details. 

The average quadratic displacement between times and 
t is approximated to the order by 



{xo[tf)^{xo[t)\(l-{sf^,,[x,p])J 



^oitfS^t,2^x,p])^, (60) 



where xo(f) - x[t) - x[0) and (•>o denotes the average over 
the action of the pure Brownian motion So [x, p] . It is shown 
in (l7j that 



(61) 



thus we only need to compute (^xoit)^ Sf^^ 2^^' Pi) The fol- 
lowing averages are needed: 



2„(f'^f,iYiQHx{t')-x[t")]^ 

X L([0, f[n[f", r'DrCiqe-^-l-'®'!*'-^"!, (62) 



and 



xo(r)V(f')p(f")Te'>'ro)-i-('''-(*")i)^ = 

Xio.tiit') [-2x[o,t[(?") + 4D;,L([0, r[n[f", r'[)r(Q)r(Q)^ 

^ p-D,riQf\t'-t"\ 



(63) 



Xi(t) is the characteristic function of the interval I cU and 
L(7) is its length. We introduce an effective damping for the 
mode Q: 



(64) 



Using these expressions and computing the integrals over 
the interaction times in the long-time limit, we get the ef- 
fective diffusion coefficient: 



In our case, it reads 



DeS 



2Dx [r(Q) -/(Q) + TriQf g(Q)Y{Q)] - rtr(g(Q))o)(Q) 



w(Q)^ 



dQ. 



h^Djc j {k+k'fk{k)^K{k'f ([2C- IJjc^ [R[k)K[k) + R[k')~Mk')] +DAk+ k'f) dfcdfc' 



2d 



A(fc)A(fc') [k^ [R{k)K{k) + R{k')Mk')] + D^[k+k']^y 



[2n) 



2d 



(65) 



(66) 



This expression, giving the effective diffusion coefficient 
of a particle quadratically coupled to a Gaussian fluctuating 
field outside the adiabatic limit is our main result; we now 
discuss some of its features. Firstly, it looks very similar to 
the expression obtained for a linear coupling in 1 17|, that is 
not surprising since the effective evolution equations, that 
are the starting point of the calculation, are very similar. The 
main qualitative difference is that the effective diffusion co- 
efficient Dx - TKx stands in front of the correction, whereas 
it is the mobility Kx for a linear coupling. This is a hallmark 
of fluctuations induced effects, where the temperature plays 
a central role. It was already observed in [18) where the drag 
coefficient was found to be proportional to the temperature. 

Secondly, the feedback parameter ( plays the same role as 
in the linear coupling case, and the consequences are thus 



the same: the diffusion coefficient of an active particle is al- 
ways reduced by its coupling to the field, whereas for a pas- 
sive particle it is reduced in a slow field and increased in a 
fast field. This behavior is represented on FIG. [2] where the 
correction to the diffusion coefficient is plotted against the 
parameters ( and k 0. 



IV. NUMERICAL SIMULATIONS 

We perform numerical simulations on a simple system, in 
dimension d-l and with a finite number of Fourier modes. 
All the parameters are kept constant except (, that is used 
to investigate active and passive diffusion, h, that allows to 
test the validity domain of the weak coupling approximation 
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Deff - D, 




FIG. 2: (Color online) Predicted effective diffusion coefficient as a 
function of the feedback parameter ( and the field evolution rate 
K^, for the parameters used in the numerical simulations. The grey 
plane indicates the value and the white line is the correspond- 
ing contour. This graph shows that the diffusion coefficient is in- 
creased for a passive particle in a fast field, and decreased other- 
wise. 



used for the predictions, and tc^, that allov^rs to vary the field 
evolution rate and thus interpolate between the quenched 
field and adiabatic limits. 

The fixed parameters and operators are set to T = 1, tCj^ = 
1, K{k) = Rik] = 1, A(fc) ^ + We use the method intro- 
duced in [21 1 to simulate a Gaussian field writh a finite num- 
ber of Fourier modes, albeit in dimension d - I here. The 
field reads 

N 

(p{y, t] ^ aoit) + ^ [at{t)cos[ky] + bkit]smiky)]. (67) 
fc=i 

To determine the evolution of the coefficients flfc(f) and 
bjcit), we insert this decomposition in the evolution equa- 
tion |8), then multiply it by cos(fcy) or sin(fcj/) and integrate 
over [0, 27i] . We get 



doit) - -K^ 
dkit) - -K(i> 

bkit) = -Jf0 



AiO)aoit)+^—(l)ixit),t) 
In 

- , C'z 

A(fc)a|fc(f) -I- — </)(x(f), t) cos(A;x(f)) 
n 



Aik)bkit] + — (pixit), f)sin(fcx(f)) 

71 



(68) 



(69) 

+ v^^'fc(f) 
(70) 



where the noises have the following correlators 



i^kitHtit'])^ 


2T ^ 
— Sit 

71 


-t')8k-k' if k^Q, 


(71) 


<^o(0^o(t')> = 


-Sit- 

71 


t'], 


(72) 


(,^',(r)4(f')> = 


2T ^ 
— Sit 

71 


- t')Sk-k', 


(73) 


i^km'tit'])^ 


0. 




(74) 
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FIG. 3: (Color online) Effective diffusion coefficient of an active 
particle as a function of the coupling parameter; simulations (cir- 
cles) and analytical result (continuous line). 



In the simulations, we let the system evolve for a long 
time T » )c~\)c|^\ and we measure the particle position at 
regular intervals. We repeat this simulation a large number 
of times (around 10^) and, using these measurements, we 
compute ( X ( f ) ^ ) , where t is the measurement time. We then 
plot this function of f and a linear fit gives the value of the 
effective diffusion constant. 

For a finite number of modes, the integral is replaced by 
a sum in the effective diffusion coefficient, and the simula- 
tions are compared to the following formula 



87r2 



N 



(75) 

.h2\ 



E 

k,k'=-N 



ik + k'f {[2( - 1] JC0 [A(fc) + A(fc')] + [k+ k'f) 
A(fc)A(fc') [k^ [\{k) + A(fc')] + [k+k']^f 



The validity domain is questioned in FIG. [s] for the active 
case and in FIG.|4]for the passive case. For an active particle, 
it appears that the validity domain is very small: the small 
coupling approximation holds for ^ < 0.7, where the correc- 
tion is only 1% of the bare diffusion coefficient. The situa- 
tion is different for a passive particle, especially in a fast field 
(jC0 - 2), where the approximation holds for h <2, where 
the correction is close to 5%. The validity domain is smaller 
for a slow field (jc^ - 0.5). We can conclude that our pre- 
diction is better when the diffusion coefficient is increased 
than when it is decreased; in the last case it rapidly under- 
estimates the diffusion coefficient. 

The effect of the field evolution rate is studied in FIG. |5] 
for a small coupling and in FIG.[6]for a high coupling. In the 
first case, the computed diffusion coefficient is close to the 
one obtained in the simtilations. We notice important vari- 
ations on simulation results; they are due to the fact that the 
correction to the bare diffusion coefficient is very small, and 
thus needs a lot of computation time to be determined pre- 
cisely. At high coupling, the agreement between theory and 
simulations is not good for an active particle, but the gen- 
eral tendency is correct. On the other hand, the passive cal- 
culation appears to be very accurate and fully reproduces 



8 



1.04 



;0 = O.5 

;0 = O.5 

= 2 
= 2 



Q 

OJ 

Q 



0.96 



0.92 



1.1 



0.5 



1.5 

h 



2.5 



FIG. 4: (Color online) Effective diffusion coefficient of an passive 
particle as a function of the coupling parameter, for two different 
field evolution rates; simulations (circles and squares) and analyt- 
ical result (lines) . 
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FIG. 6: (Color online) Effective diffusion coefficient as a function 
of the field evolution rate for a passive and an active particles, with 
a coupling h = 2; simulations (circles and squares) and analytical 
result (lines). 
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FIG. 5: (Color online) Effective diffusion coefficient as a function 
of the field evolution rate for a passive and an active particles, with 
a coupling h = 0.5; simulations (circles and squares) and analytical 
result (lines). 



the simulations results. Notably, the crossover between a 
decrease and an augmentation of the diffusion coefficient 
when the field becomes faster is well captured by the ana- 
lytic calculation. 

As a conclusion, our predictions are precise in a relatively 
small domain of validity, but they remain qualitatively good 
at higher coupling, and thus provide a useful understanding 
of the phenomenon studied here. 



V. DISCUSSION ON ANOMALOUS DIFFUSION 

The study presented here focuses on the computation of 
small corrections of the effective diffusion constant in the 
case where the diffusion remains normal. However, as ex- 
plained in (ITj 1^ an infrared (small k) divergence in the 
integral invoked in the expression of the diffusion constant 



indicates a fertile ground for anomalous diffusion. We now 
discuss the infrared divergence of the integral appearing in 
66 , in a way similar to the one used in [17|: we start by 
defining the small k behavior of the operators. 



A(fc) 



k', 



ifc-O 

K{k] ~ k", 
fc— 

R{k] ~ kP. 



The integral in 66 diverges at low k if 

dego = 4k - 25 - min(p -1- 5, 2) -1- 2d -1- 2 < 0, 



(76) 
(77) 
(78) 

(79) 



where dego '^he degree of the integral, i.e. its dimension 
in k. In the two examples presented in section [ll] the op- 
erator A is of the form A[k) - k^ [k^ + m^). When the field 
is critical, m-0 and S' - 6 + 2 should be used instead of S 
in the previous formula. We now apply this formula to our 
examples. 

The following table gives, for the cases of a protein cou- 
pled to membrane curvature or composition, the value of 



the exponents and the degree dego of the integral 66 





curvature 


composition 


parameters 


5 = 2, k: = 2, p = -l 


5 = 0,)c = 0, p = 2 


non critical 


dego = 9 


dego = 4 


critical 


dego = 4 


dego = 



It appears that a coupling to membrane composition at 
the miscibility transition can lead to anomalous diffusion. 
Since the proteins are active {( - 1), the correction to the 
diffusion constant is negative and the anomalous diffusion 
should be subdiffusive. We recall that this statement should 
be taken carefully: we only say that this system is likely to 
be subdiffusive, we do not assert that it is. For instance, this 
divergence could also stand for a correction proportional 
to h°, with a < 2. The degree of the integral implies a 
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logarithmic divergence, i.e. a correction of the form Iog(/) 
when the system is close to criticality. In this configuration, 
the temperature dependence of the diffusion coefficient is 
dominated by the temperature dependence of the correla- 
tion length ffg) . 

We now rapidly discuss the case of ultraviolet (large k) di- 
vergence. An ultraviolet divergence should be regularized 
with a microscopic cut-off length, for instance the protein 
size; in this case the diffusion coefficient thus depends on 
the protein radius, an effect that was explored experimen- 
tally in |6| and discussed in [14l[T8]. The divergence crite- 
rion is here 

deg^ = 4?c - 25 - max(p + 5,2) + 2d + 2>0, (80) 

where the exponents are defined in the limit k (x>; this 
analysis does not depend on the correlation length. We find 
that deg^ - 3 for a coupling to the curvature and deg^^ - -2 
for a coupling to the composition. In the first case, regular- 
izing the divergence with the protein radius a gives a cor- 
rection proportional to a^. 

VI. CONCLUSION 

We analyzed the diffusion of a particle quadratically cou- 
pled to a fluctuating field, that avoids regions where the field 
fluctuates and may reduce the field fluctuations, in which 
case it is called active; if it does not affect the field it is 
called passive. This study completes the one developed in 
17 1 for a linear coupling between the field and the parti- 
cle and they set together a general model for diffusion in a 
fluctuating environment. The environment is modeled by 
a Gaussian field obeying overdamped Langevin dynamics. 
The particle-field interaction can take two forms: if the par- 
ticle breaks the symmetry of the field and prefers one spe- 
cific sign, it is modeled by a linear interaction; if it does 
not break the symmetry of the field and avoids field fluc- 
tuations, it is represented by a quadratic interaction. Both 
types are present in the general case, an example being the 
protein coupled to membrane curvature discussed in sec- 
tion |IIB| In the limit of a weak interaction, both interaction 
types contribute to the full correction to the effective diffu- 
sion coefficient; alternatively a preliminary study may de- 
termine the dominant contribution. 

We focused on the computation of the effective diffusion 
coefficient, assuming normal diffusion; our method allows 
to deal with equilibrium as well as out of equilibrium sys- 
tems (that is the case if the environment is not affected 
by the particle or subject to non thermal forces |3)|29|). 
Our conclusion is that for both types of coupling, the dif- 
fusion coefficient is always reduced for an active particle, 
whereas it can be reduced or increased for a passive particle 
if the field is slow or fast, respectively; this is summarized in 
FIG.H] 

When the correction to the bare diffusion coefficient di- 
verges, it may be interpreted as the onset of anomalous dif- 
fusion; depending on the sign of the correction, it indicates 
subdiffusion or superdiffusion. For instance, for membrane 



proteins, we have shown that subdiffusion may occur when 
the protein is coupled to the membrane composition at the 
critical point. This output of our theory is of major im- 
portance regarding the numerous examples of anomalous 
diffusion found in experiments |2, 3); however, it gives lit- 
tle information on the diffusion and notably on the expo- 
nent. This limitation comes from the fact that we use a per- 
turbative expansion around normal diffusion and compute 
the prefactor, i.e. the diffusion coefficient. Further work is 
thus needed to be able to deal with environment induced 
anomalous diffusion; a necessary step in this direction is to 
go beyond perturbative results. 
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Appendix A: Functional operators in real and Fourier space 

In this appendix, we define our notations and recall some 
basic properties of functional operators. We start in real 
space, and then see how it transposes to Fourier space, that 
is used a lot in this article. All the operators considered here 
are real. 

For two functions /(y) and g{y) and two operators 
A{y,y'] and B{y,y'], the scalar product of / and g, the ac- 
tion of A on /, the product of A and B and the trace of the 
operator A are respectively defined by 





1 fiy)giy)dy, 


(Al) 


(A/){y) = ^ 


f A{y,y')f{y')dy', 


(A2) 


(AB)(j/,y')=^ 


f A{y,y")B{y",y')dy", 


(A3) 


tr(^) = 


/ A{y,y)dy. 


(A4) 



The adjoint A^ of the operator A is defined by the equality 
/ • [Ag) - [A^f] ■ g, that holds for all functions / and g; it is 
straightforward to show that 

Ahy,y')^A(y',y). (A5) 

A is symmetric or self-adjoint if A^ - A, or 

A[y,y')^A[y',y). (A6) 

An operator A is invariant by translation if there exists a 
function a such that 

A{y,y')^a[y-y'). (A7) 

Such an operator is isotropic if it only depends on the dis- 
tance between y and y': 

A[y,y')^a{\y-y'\). (A8) 
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We now switch to the Fourier space, with the Fourier 
transform defined by 



dk 



dfcdfc' 



(A9) 
(AlO) 



for a function and an operator, respectively. This definition 
allows us to translate |A1[|A4) into Fourier space: 



fi-k)g{k) 



dk 



Af[k)^ j A{k-k')f[k')^^, 



AB[k,k')^ \ A[k,-k")B(k" ,k') 



dk" 



{2n) 



d ' 



ti{A)^ A{k,-k)- 



dk 



(All) 
(A12) 
(A13) 
(A14) 



The adjoint of the operator A is defined in the same way as 
in real space: 



AHk, k')^A[k',k). 



(A15) 



The Fourier transform of the translation-invariant operator 
My,y') - a{y- y') reads 



A{k,k') ^ {Znr dik)5{k+ k']. 



(A16) 



Moreover, if A is isotropic, its Fourier transform only de- 
pends on the norm \ k\: a{k) - a{\k\). 

In this article, we do not use a different notation for the 
one-variable function associated to a translation-invariant 
operator: the number of variables indicates if we refer to the 
operator or to its associated function. For instance, we will 
use A(fc,fc') = (.27i)'^Aik)S{k+k']. 



Appendix B: Field correlation functions 

In this section, we derive the free field correlation func 
tions (y/if), (cpy/) and (cfxp) from the free field action 22 



Equations on averages can be obtained from the action with 
the Schwinger-Dyson formula 1 30 1 : 



C-^)JAicf>,^,<^ 
\ diPiy, t) I \ ^ ^ dcpiy, t) 



(Bl) 



where A is a functional operator. The same equation holds 
if the derivative is taken with respect to i//(y, t). Using ^ = 1, 
we derive two equations that give the averages 



{ip(y, t)) = 0, 
<i/'(y, f)> = 0. 



(B2) 
(B3) 



With A[(p,\j/] - y/{y,t) and derivating with respect to 
(p{y', t'), we get 

= {y/[y, t) l^friy', t')-Ky,{RAy/)[y' , t')]), (B4) 
that leads to 

C„iy,y', t, t'] = {yj[y, t)^p{y' , f')> = 0. (B5) 

Setting A[^,%ir] - (p{y,t) and derivating with respect to 
^iy' , t') gives 

8{y-y')8(t- t') = i{(p{y, t) [y;(y' , t')-K^iRAy/){y', t')]). 

(B6) 

This equation has the solution 

C^y,iy,y', t, t') = (cpiy, t)rfi[y', t')) (B7) 

Finally, keeping the same operator A and derivating with re- 
spect to V(y'. f')' we get 

i((/)(y, t) [<p(y', t') + K^[RA(P)(y', r')]) = 

2TK^{(p{y,t)(Ryj)[y',t')). (B9) 

Using the previous result ]B8) , we solve this equation with 



C00(y,y,f, f') = (0(y, m{y',t']) 



(BIO) 

= r[A-ie-'^*l*-^'l^^](y-y'). (Bll) 
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